Depinning by Fracture in a Glassy Background 
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We force a single particle through a two-dimensional simulated glass. We find that the particle 
velocity obeys a robust power law that persists to drives well above threshold. As the single driven 
particle moves, it induces cooperative distortions in the surrounding medium. We show theoretically 
that a fracture model for these distortions produces power law behavior, and discuss implications 
for experimental probes of soft matter systems. 



PACS: 82.70.-y,74.60.Ge,64.60.Ht 

Recently, there has been considerable interest in the 
behavior of systems of interacting particles that resist 
flow, including vortex glasses in superconductors [1,2], 
charge-density- wave metals [3] , Wigner crystals [4] , dense 
colloidal suspensions [5,6], and Coulomb-blockade arrays 
[7,8]. In each of these systems, quenched disorder and 
repulsive interactions from the surrounding particles pre- 
vent any individual particle from moving in response to a 
small external force. As a result, below a threshold force, 
the system behaves as a solid, while above the thresh- 
old, the system can flow plastically and is soft and disor- 
dered. Once the system depins, increasing the force leads 
to rapidly increasing motion, as the speed and transport 
response increase faster than linearly with the flow over a 
range of applied force. Theories [7,9-11] have attributed 
this anomalous transport to the concerted action of many 
driven particles across a rough pinscape. They predict 
power-law growth of the velocity with force in the vicin- 
ity of a threshold force Fc- v ~ {F—Fc)^. In experiments 
and simulations, the velocity above threshold has been 
observed to vary as a power /3 of the force, with (3 in the 
range 1.5 to 2.2 [4,5,12,13]. In contrast to the critical- 
state model of elastic depinning [9] , where scaling occurs 
only very close to threshold, the power law is observed 
to hold for forces of several times the threshold force. 

In this work we demonstrate power-law collective 
transport with a single driven particle in a disordered 
glassy matrix of other non-driven particles in two dimen- 
sions, realized via a molecular-dynamics simulation. We 
flnd a power law with /? = 1.5 over two decades of force 
that is insensitive to the system size, the size of the driven 
particle, or the density of the surrounding medium. This 
contrasts with a single particle driven over a substrate 
with quenched disorder where a scaling of /? = 1/2 is 
expected [9] . Our result suggests that the origin of the 
(3 > 1 scaling in a variety of systems may be simpler 
than previously supposed. A single driven particle drags 
other particles with it, thus slowing it down. The faster 
it moves, the fewer particles it drags, and hence the expo- 
nent P becomes larger than 1.0. We explicitly show that 
an analysis in terms of fracture in front of the particle 
gives (3 = 1.5. The fracture leads to a one dimensional 
(ID) plastic zone, which appears as a river-like flow of 



particles [2]. While in other cases such rivers are at- 
tributed to easy paths through a background of quenched 
disorder, in our system the only disorder is due to the 
glassiness of the medium. The plastic zone appears due 
to the softness of the system; for stiffer interactions, the 
plastic zone disappears and P ~ 1 is the observed scaling. 

In addition to offering a simpler model of depinning, 
our results are relevant to systems in which two species 
of particles move in opposite directions with respect to 
each other. Such situations can arise in certain types 
of electrophoresis experiments, pedestrian motion, self 
driven particles and molecular motors [14], as well as in 
vortex experiments in which a vortex lattice is driven 
past a low density of pinning sites. 

Further, driving a single particle through a soft matter 
system can be used as a powerful experimental probe of 
dynamics of the medium far from equilibrium. For exam- 
ple, in recent experiments a magnetic particle is dragged 
through a colloidal system near the glass transition [6]. 
The absence of momentum conservation in our system 
leads to very different physics from the transfer of mo- 
mentum through hydrodynamic and hard-core interac- 
tions. However, such experiments on colloidal systems 
and emulsions in conflned geometries, in particular be- 
tween parallel walls, are predicted to have similar fea- 
tures to those found here. 

Simulations and Dissipation Balance- We drive a sin- 
gle particle with a charge at constant force through 
a two dimensional (2D) disordered system with periodic 
boundary conditions in the x and y directions. To create 
a glassy medium and prevent formation of a triangular 
lattice, we simulate a mixture of two species of particles 
with different charges qa = q and = 2g, where g = 1. 
For the driven particle, we consider a range of values from 
qo = 0.15g to qn = BOg. The overdamped equation of 
motion for particle i is 77V = f; = — J2j ^U{rij) + F^x, 
where v is the particle velocity, the damping coefficient 
r] — 1, and = on all particles except the one with 
charge qo- We use a screened Coulomb interaction, given 
by: 



Uirij) = q,qj- 



-2rij 



(1) 



1 



between particles i and j separated by a distance Vij. 
We have considered a variety of system sizes for N = 480 
and N = 2150 particles, as well as lattice constants of 
a = 1.08 to 2.3. 

In Fig. 1 we show the positions of the particles and 
the trajectories for a fixed number of time steps for a 
system with = 2100 particles at a drive Fj = 0.75. 
The driven particle, marked as a large dot, has qn = 20. 
We find that in general there is a finite threshold force 
Fc for the particle to move. The perturbation of the 
other particles by the driven particle is anisotropic, with 
a larger perturbation in the direction of drive than in the 
transverse direction. Particles more than a few a away 
from the path of the driven particle move elastically in 
small nearly closed orbits of radius less than a. Particles 
in front and behind the driven particle exhibit plastic 
motion. This plastic zone decreases in size for higher 
velocities. 

The medium, while disordered, is a solid. It has a 
threshold, F^ = Fc, for failure, and can support shear 
stress. Thus, the particle must fracture the medium to 
move through it. From momentum balance, we have that 
at all times F^ = rjV + 'qY^VA, where V is the velocity of 
the driven particle and the sum ranges over all other par- 
ticles. Below threshold, where the medium moves with 
the particle, the momentum balance yields V = F^/N, 
where A'' is the number of particles. Above threshold, to 
obtain power law scaling, some large number n of other 
particles, 1 <C n <C VV, must move with the driven par- 
ticle, with some finite size corrections to scaling due to 
net background motion of the remaining N — n particles. 
Eventually, as n — > 1 far above threshold, the velocity 
returns to linear scaling. 

In Fig. 2 (lower solid curve) we show the log-log plot 
for the velocity vs applied drive F^ — Fc, for the driven 
particle in Fig. 1, where Fc = 0.7 is the threshold for the 
large particle to move. We find a good scaling with a fit 
of /3 = 1.47 ± 0.03. We find a similar scaling for a variety 
of other parameters with the exponents of 1.5 ± 0.05. At 
sufficiently large drive, /3 returns to 1, as expected. The 
size of the scaling region decreases with decreasing qo /q, 
until for qo/q < 1-0, no anomalous scaling is observed. 
In this case the driven particle does not induce plastic- 
ity in the other particles but only a smaller perturbation 
of size less than a. In the upper solid curve in Fig. 2 we 
show the velocity force curve for = 0.5 showing a scal- 
ing fit with (3 = 1.0. For a fixed value of Fd the velocity 
decreases with increasing q^ as the driving particle in- 
teracts more strongly with the siirrounding particles. In 
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Fig. 3(b), V vs qD for F^ = 4.0 shows a g^, scaling in 
the regime where there is plastic deformation, and then 
V flattens out in the elastic regime ior qo < 1. Addi- 
tionally in the plastic flow regime the particle motion is 
highly intermittent as illustrated in Fig. 3(a), where a 
time trace of V for qn/q = '20.0 shows the motion occur- 



ring in bursts separated by quiet periods. As the drive 
increases, the motion becomes less intermittent, and the 
noise spectrum of the time trace develops a well-defined 
rollover frequency w, with cj scaling approximately lin- 
early in Fd — Fc- We find that the rollover frequency 
corresponds to the time scale between peaks in Fig. 3(a), 
rather than to the duration of a single peak. 

In Fig. 4, we show a contour plot of the average of v"^ 
for the other particles in the medium, as a function of 
position relative to the driven particle, for systems with 
qo/q = 20, N = 2100, a = 1.085 and Fd = 0.9,1.5, 
and 40 respectively. This measures energy dissipation, 
as from energy balance FdV = rj{V'^ + X^'^a)- These 
systems are above the depinning threshold of 0.7, but 
still within the scaling regime. The dissipation is cen- 
tered strongly around the moving particle, and extends 
anisotropically into the surrounding medium, with a 
larger region of dissipation in the direction of drive. 

Theory: Elasticity and Fluid Flow We now consider 
various scenarios for the particle motion to understand 
the simulation results. We begin by considering the re- 
sponse of the medium away from the plastic zone, where 
the particles move less than a and elasticity theory is ap- 
plicable. The driven particle exerts a force on the elastic 
medium at a location that is moving at velocity V . At 
long wavelengths, elasticity theory gives for a lattice dis- 
placement u(k) the equation of motion 

U{k) = Cik{k ■ u) + C'ik X {ux fc). (2) 

To obtain elastic constants Ci, C2, we consider a triangu- 
lar lattice of particles, with appropriate lattice constant 
a, and assume that all particles have the same charge, 
qA,B = V2q. For a lattice constant a = 1.085 and the 
chosen q w 0.086, one finds Ci w 0.025/r7, C2 « 0.0018/7?. 
The particle exerts a force on this medium, along the di- 
rection of motion of the particle. There is also a force 
normal to this direction, pushing the medium out side- 
ways in opposite directions on opposite sites of the mov- 
ing fracture. In the co-moving frame, the displacement 
in any given direction decays exponentially at large dis- 
tance with a characteristic length I = Ci/V along that 
direction and I = C2/V normal to the direction. The 
energy dissipation rate resulting is of order rjV (l/a)^ = 
r]CiC2/a^, and is thus independent of V. This indicates 
a force on the driven particle of order 1/ V, in addition to 
the drag force, r]V , on the driven particle. For small V, 
this force becomes arbitrarily large, ultimately exceeding 
the yield stress of the medium. In such a region the ma- 
terial must yield and the elastic picture becomes invalid. 

Thus, we must consider the plastic zone. It is not 
consistent to have a 2D chunk of solid moving with the 
particle, with a plastic zone between that solid and the 
rest of the medium: if the particle force is sufficient to 
fracture the medium, it would fracture the solid moving 
with it. At the same time, one cannot have a 2D fluid 
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zone around the particle. The overdamped nature of the 
dynamics (assuming viscosity and pressure terms hke a 
Newtonian fluid) would concentrate all the vorticity of 
the fluid within a length a around the particle, and thus 
this scenario is also inconsistent: more than a length a 
from the particle, the stress is too small to fracture the 
medium and produce the fluid. By elimination we turn 
next to a ID plastic zone. 

Theory: Compressed Column — We now propose a sce- 
nario based on competition between shear and compres- 
sive failure which accounts for all the numerical results, 
and illustrate it using the specific systems of Fig. 4. The 
driven particle creates a fracture in front of it. Naively 
the fracture force would be expected to be of order the 
particle interaction force aU"{a), or roughly 0.015 for 
the given qd = 20. In fact, it is equal to 0.7, nearly two 
orders of magnitude greater, for two reasons: the force 
U'{r) rises rapidly at short r, while the charge Qd > Q 
further increases the fracture force. After fracture, the 
particles in front of the driven particle must then move 
out of its way, either by failing in shear and moving along 
with the driven particle, or by failing in compression and 
moving transversely out of its way. Initially, consider just 
the first possibility so that in front of the driven particle 
there is a growing ID column of particles failing in shear. 
One finds that the rate at which this length increases 
is determined not by the velocity of the driven particle, 
but by the (faster) velocity at which a compressive front 
ahead of that particle moves. ID elasticity theory would 
imply that in time t, this compression reaches a distance 
I « a^JVi/rj. Since this system is far from equilibrium, 
we have checked this result by simulating a ID system 
with a single driven particle in a fixed, periodic back- 
ground potential to mimic the medium. We find that 
this behavior remains valid, even close to the depinning 
transition, albeit with a greatly increased value for U" . 

The ID column does not grow indefinitely in length due 
to the possibility of transverse motion of the particles, in 
which case they squirt out of the column. The time scale 
for this process would naturally be of order rja/ (Fd — Fc). 
Thus, the number of particles in the column is l/a « 
\/U"a/{Fd — Fc). By momentum balance, this number 
is equal to VjV /{F^ — Fc), giving 

vex (i^d-F,)3/2, (3) 

as observed, where some fixed force Fc is required to cre- 
ate the fracture. The exit of particles from the column 
due to compressive failure is naturally highly intermit- 
tent, leading to intermittent motion of the driven parti- 
cle, as seen in Fig. 3(a). This predicts the hnear scahng 
of the rollover frequency with F^ — Fc, with the prefac- 
tors such that the rollover frequency is increased above 
the naive expectation: at F^ = 1.5, ui = 0.05, larger than 
rja/iFd - Fc). 

The ID plastic zone is visible in the contours shown 
in Fig. 4. The aspect ratio of the contour increases for 



lower drives as the length of the contours increases at 
constant width, confirming that dissipation arises in a 
ID plastic zone, rather than due to the elastic response 
of a 2D medium which would instead give a constant 
aspect ratio. Further, the size of the contours is not 
consistent with 2D elastic response, which would appear 
in Fig. 4(b) as an exponential decay of the dissipation on 
a length scale of approximately 0.35 for V = 0.072, while 
the actual contours are significantly larger. 

As the charge of the driven particle increases, it pro- 
duces a depletion zone of missing particles around it to 
compensate its charge; the area of this zone is propor- 
tional to qjj I q. The length of the compression was deter- 
mined above; the width will be proportional to the radius 
of the depletion zone, ^/qo/l-, accounting for the scaling 
of velocity with q^, as seen in Fig. 3(b). 

Summary — We have found a robust power-law for the 
velocity of a single driven particle fracturing a glassy en- 
vironment, with /? = 1.5 ± 0.05. We give a theoretical 
explanation based on a ID plastic zone, which decreases 
in size as the particle moves faster leading to /3 > 1. This 
behavior arises due to the softness of the material, as for 
larger lattice constants a or smaller qu, the interaction 
becomes stiff'er, and the width of the scaling region de- 
creases until it eventually disappears leaving only linear 
scaling. While for conventional solids this scaling region 
is absent with no column observed in front of the particle, 
the origin of anomalous transport features in many disor- 
dered systems may be due to the mechanisms discussed 
herein. 
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FIG. 1. The particle positions (black dots) and trajecto- 
ries (black lines). The driven particle (large dot) has charge 
qo = 20; the remaining particles are a mixture oi qA = q and 
qB ~ 2q. There are a total of 2100 particles. The trajectories 
are drawn for a fixed number of time steps with a constant 
applied drive of Fd = 0.75. 
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FIG. 2. The velocity V vs Fd — Fc for (upper solid curve) 
a driven particle with qo = 0.5. The upper dashed curve is a 
fit with (3 = 1.0. The lower solid curve is for qo ~ 20, with 
the lower dashed curve a fit with /3 = 1.47. 
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FIG. 3. (a) The velocity vs time for a fixed drive of 
Fn =0.75 with qo =20, showing intermittent bursts of mo- 
tion, (b) The velocity vs qn for a fixed drive Fd = 4.0. The 
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solid line is a fit of . 
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FIG. 4. Contour plot of as a function of position relative 
to the driven particle. The lowest contour is normalized to 
5% of the maximum of , and each contour represents a 10% 
increase. The system size is 48 x 48, and the driven particle 
position is fixed to (24, 24). 
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